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Introduction 


Calculation of transition matrix elements between states represented by Cou- 
lomb wave functions are important in many physical applications. Special cases 
of such matrix elements were considered by W. Gordon [1] and more recently 
K. Alder and A. Winther [2] have given general expressions and recursion for- 
mulas for the matrix elements, which are given in terms of hypergeometric 
functions of two variables. These functions are in many cases comparatively 
difficult to compute due to the lack of representations of the functions that 
lend themselves readily to computations Their series expansions often converge 
for limited values of the occurring parameters only and this restriction is some- 
times a severe one. Integral representations often have the complexity of singular 
integrands unless detours of the path of integration are made into the complex 
plane. When relativistic effects are of importance, the non-integer steps of cer- 
tain parameters decrease the usefulness of recursion relations and, as cross-sections 
may depend on a large number of matrix elements, it may be necessary to 
compute a considerable number of them directly. 

In a previous paper by the present author [3] some aspects on Coulomb wave 
functions were given. The wave functions were expressed in terms of generalized 
hypergeometric functions in an attempt to increase the number of representations 
of the matrix elements as a result of the introduction of more general functions. 
As will be shown in this paper, the expressions obtained in [3] lead rather di- 
rectly to simple expansions of certain matrix elements that are convergent for 
all values of the occurring parameters that have a physical meaning. The con- 
vergence of the expansions that will be presented is fastest for small energies. 
tendard treatment of the considered matrix elements leads to expansions that, 
for practical purposes, are best convergent for high energies. 


Transition matrix elements between bound and non-bound states 


‘Higher order perturbation calculations involve matrix elements between arbi- 
Coulomb states or, more precisely, the calculation of matrix elements for 


159 


PER 0. M. OLSSON, On transition matrix elements between states in a Coulomb field 


transitions between non-bound states cannot be avoided. These matrix elements 
are sufficiently complicated to make such an undertaking as e.g. a rigorous 
calculation of Bremsstrahlung a very tedious matter even if a fast computor is 
used. The matrix elements entering into a first order pertubation calculation 
for transition from a bound to a non-bound state are simpler, particularly if 
the interaction may be written essentially as an exponential function as is the 
case when radiative interaction is caused by photons having well defined angular 
momentum. An example is internal conversion. In the similar process of ex- 
ternal conversion or the photo effect, the incident photons must be described 
by a plane wave. The spherical symmetry of the problem necessitates expan- 
sion in multipoles if partial waves are to be used. The matrix elements en- 
tering into calculations of photo cross-sections are then of the same type as 
the matrix elements entering into calculations of internal conversion coefficients 
but they are more difficult to compute due to the fact that there is no pre- 
dominant near-zone of the radiation field as is the case in internal conversion, 
where the radiation is emitted from the nucleus. 

The following considerations will apply to matrix elements consisting essen- 
tially of integrals of the type: 


eo 
I= I et 7t-t otkr g-ier BL (y+in; 2y; 2ipr)dr. (1) 
0 
Here e “’ and e*’ may be thought of as the exponential functions entering into 
a bound wave function and a radiation field respectively. All parameters are 
intended to be real and positive. The function 


-ipr 


eo" Fi (ytin; 2y; 2tpr) (2) 
is real valued and stands for a non-bound wave function, whereas r*? contains 
the occurring powers of r. Real expressions for the non-bound Coulomb wave 
functions, which are usually given as real and imaginary parts of a complex 
valued function, have been given in a previous paper [3]. The integral-(1) is 
a special case of a more general integral computed by A. Erdélyi [4]. Term- 
wise integration yields 


a i e Meter el BP (y+in; 2y; 2ipr)dr= 


= $ (y+tn)n coEcE eo (Ho ik+ipyr ptt L Ties 
n=0 (2y)nn 
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n=-0 = (2) nM (u—th+ip)e*” 
mony 2ip 
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We have here used the power series expansions 


; , = (ytty)n (2tpr)” 
. . 2 ay 
Ff, (y+in; 2y; 2ipr) 2, em 


a siete 2Qip @ 
ee (2y)n 7 SAS eee 2) 


(4) 


iy 2ip 
F(a yin 2 P| = zi 


and the symbols 
(y tin) =(y+in) (y+1+in) (y+2+in) ... (y+n—1+in) ete. 


I'(q) is the [-function of argument qg. The series expansion (5) for the hyper- 
geometric function in (3) converges when 


ren a 
|u—iktip|~ eo 


The fact that the function (2) is real valued implies that 
e'" Fy (y tin; 2y; 2ipr) =e” Fy (y—in; 2y; —2ipr) 


and thus the integral (1) does not change its value when the signs of p and 7 
are changed simultaneously. 
Then 


[ (q) +2ip ). 
_ + 
oe ae oa MORES Sieg ae 


If the power series expansion is to be used for computing the matrix ele- 
ment, it is natural to choose the sign that makes the modulus of the argument 
smallest. Denoting the argument by z we have 


. = 
ee A 
in p—t(k+p) (8) 


\ 
} 


Let us investigate to some extent the possibilities of the expression (7) from 
numerical point of view. A reasonably fast convergent power series expansion 
beyond doubt the most convenient way of computing it except for special 
alues of the parameters. It is convenient to have a special problem in mind 
and let us as an example choose internal conversion. Then the integral (7) is 
rt of the expression for the probability that a bound electron is converted 
a free electron with momentum p as a result of a nuclear transition. A 
peting process is the emission of y-rays with wave number k. For con- 
ion in the atomic K-shell we have 4=«Z, where « is the fine structure 
mstant and Z is the nuclear charge number. The other parameters are 
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y=V1-# 2 
_HaZ 
4 Pp 


E=V)1+>". 


We are here considering relativistic wave functions. The parameter q depends 
on the nuclear spin change but will have no influence on the following argu- 
ments. Since y-ray emission and conversion are competing processes following 
the same nuclear transitions we have 


k=H—-y 


expressing that the y-ray energy is equal to the kinetic energy of the conver- 
sion electron plus its binding energy. Figure 1 shows the argument 2 as a func- 
tion of k in the complex z-plane for Z=82. The trend of the curve is the 
same for all Z values. 

The figures along the curve give the k-values in units of mc? =511 key. 
It is seen that the power series expansion (4) of (7) converges for k<0.3 ap- 
proximately where the curve cuts the unit circle. This corresponds to y-ray 
energies less than 150 kev and conversion electron energies less than 50 kev 
approximately. For lower Z values, the k-values for which convergence ceases 
are smaller and most nuclear transitions of interest fall beyond the limit of 
convergence. Then analytic continuations of the ,/, function in terms of other 
/,-functions can be used. E. M. Rose and collaborators used expansions in 1/z 
in their earlier calculations of internal conversion coefficients [5]. Figure 2 shows 
the reciprocal argument 


3 ce th— PB) 
2ip 

as a function of k and the expansion is now seen to be convergent for k >0.25 
approximately. : 

The fact, that z+1 when k->co in figure 1, suggests expansion in powers 
of 1—z instead of 1/2: 

2ip aLZ—t(k—p) 
Z—i(k+p) aZ—i(k+p) 


L—2=1 
a 


It is easily seen that such an expansion is convergent for all k-values except 
possibly for the threshold value for conversion, which is the point 0 in fig. 1. 
We therefore use the well known analytic continuation of the hypergeometric 
function that changes the argument from z to 1—z: 


Ue Deh (cae) 
ally (4, 0: 65 @) = Tee ay Te 0) * 


T' (c) D(a+b-—c) 
TP (a) 1 (6) 


F(a; 6; a+6—c+1; 0—2) 


+ 


(l—z)°-*~? ,F, (c—a, c—b; c—a—b+1; 1-2). 
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Fig. 1. Fig. 2. 


Then we obtain for the integral (7) 


I T(q) T(2y) | LT (y—q+tn) 
(aZ—tk—ip)* \P (2y—q)T (y +i) 


: : aZ—ik+ip 
yy at gga = 5 eae ER 
4 (a7 nn q-¥ int eae) + 


I (q-—y-—in) (Zaseeeey™. 


Dg) VT (y—in) \aZ—ik—-ip 
a 4, aZ—ikt+ip\) 
oF, (2y—q, y+in; y tin; 27a - 


The expression is relatively complicated considering the simplicity of the in- 
tegral to be computed. The expansion in powers of 1/z leads to a very similar 
expression. It seems questionable whether these expressions are more practical 
to use than the calculation of the integral numerically without using analytical ex- 
pressions. 

In a previous paper, Coulomb wave functions were expressed with the aid 
of a confluent Appell function [3]. Calculations were carried out for the rela- 
tivistic case. It was shown that the well known equations for the radial part 
of the Dirac equation: 


Se cee (« 1422) g 


dr r 

d Z 1 (70) 
g a4 Be oh 

s2_(w+1+%2) 7 re 


a - © 


have solutions 
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OZ aes Eu+y ¢ 
= lo , £,t)+ i) +1, 2, 
. wee : paket! : , (11) 
:: Ley ot 
y-e lo , #,t)+ ® +1, 8,9) 
9 Et aa Oe ) 
where t=2aZr and 
®D (y, LE, t)h=O,(y+in, y—in; 2y;tpr, —ipr). (12) 


Two independent solutions are obtained for two possible choices 
y= + V2 = oF Z. 


The function ®,(y+%in, y—tn, 2y, ipr, —ipr) is one of the confluent Appell 
functions which is defined as follows: 


(5) m (Bn 2 y” 
0 (C)min m! n! 


®, (b, b’; ¢; a, y) = (13) 


11M8 
o 


This series is convergent for all finite values of x and y. In [3] it was proved 
that the function (12) is reduced to more elementary functions when b+)’ =c 


®, (6, b'; b+.0'; x, y) =e" F, (6, 6+b'; x—y). (14) 


This connects the solution (11) with more familiar expressions. The reduction 
(14) of the confluent Appell function seems to have first been given by Burch- 
nall and Chaundy [6]. From (14), it follows easily that 


e ®, (6, b'; b+6'; x, y) =, (0, 6’; b+b'; +A, y+A), (15) 


valid for any A. Thus an exponential function as factor results in the exponent 
being added to the argument of the function. From (11) it is clear that the 
integral (1) can now be written 


Co 


I= [estate O,(y+in, y—in; 2y; ipr, —ipr)dr (16) 


0 


from which a simple expression for the integral may be derived. Using (15) 
we obtain 


Te | er D (ytin, y—in; 2y;i(k+p)r, i(k—p)r)dr. (17) 
0 


The expansion (13) is inserted and the resulting series is integrated term by 
term, which yields 
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= + inn (yim) Li (P+ DT" =D)" DO) @men 
m,n (2) men m! n! yagi 


LOM y Omen (vy tty)m(y—t)n ( eee)" (' aa 


ie m,n (2 V)m+n m! n! Lt lu 


This double series is the Appell function F, and thus 


i(k+ p) =P). (18) 


C@ 
I= TTR OS Nf og) Rae RO IIR, ed 
Feed as ‘ial hd cad AMAL Fava! r i 


The series expansion of F, used above is convergent when the moduli of both 
arguments are less than unity [7], which means 


The function F, allows, however, a transformation 


F(a [ees a ear eater ee a (ei eens UE 
14@,.0,.0; ¢; z, y)=(1—z) (1-9) P,(c CATS Sesrans rps y (19) 


and as the arguments in (18) are purely imaginary as a consequence of the 
use of the transformation (15), it is seen that the transformation (19) leads to 
an expression for the integral where the hypergeometric series defining F, is 
convergent for all values of the parameters that have a physical meaning. In 


addition, it is an expansion where the convergence is fastest for small energies. 
Using (19) we finally obtain 


on I (q) ( —1 ete" (oHe= a . 
ue lu u 


i(k+p)  _i(k—p) ) 
By (2y— ay bin yim 2; 17 - tO} 
{27% 79D YO AY Tee) — wi (k—p)— am) 
_ There exist many expansions of F', in simple series but care must be taken 
not to loose the convergence when such an expansion is used. Among the pos- 
ible expansions we select 


7 ide ia aa = (4)n a em UR {Se meg, n ¢ 
FG, bi 5'; ¢; 2, a a es a 2F,( n, b';b+6'; 1 4) ; (21) 


The derivation of (21) given by Appell and Kampé de Feériet [7] is based on 
e integral representation 
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F, (a, 6, 6’; ¢; x, y) =const. 


uae | 

| free dap ht [Lo tty tay dod (22) 
0 0 

where the last factor in the intergrand is expanded: 


{[l—s(a+ty—ta)] °= > De gn (otty—ta)”, (23) 
n=0 ‘ 


This expansion is inserted in (22). After term-wise integration, the expansion (21) 
is obtained. The expansion (23) is uniformly convergent when 


|at+ty—ta|<1, O0<t<\l. (24) 
With the arguments of the function in (20) we have 


i(kt+p) _(k+p)*—iw(e+p)_1 1 (k+ p)—p?—2t pu (k+p) 


il 

vi (k+ DP) — be e+ (b+ p) 2° 2 pe+(k+ p)? 
1 
2 


l(k+p—ip)? 1,1 
2 (k+pr+we 2 2 


elf 


where 91 =2 aretg ; - 
eee 
because (kK+p mh ae 
(k+-p)" + u 
i al 
Similarly Sees a = e*, g,=2 arctg k = p 


The two points x and y are thus always found on a circle in the complex 
plane with radius 1/2 and its center at the point 1 /2. No points on the pe- 
riphery of this circle are situated outside the unit circle and the point ~+ty—tx 
is never outside the smaller circle with radius 1/2 as it is always found on a 
straight line having x and y as end points. This proves the convergence of the 
expansion (21). 

The expansion of our expression (18) is then finally 


Fah le ectete (le 
lps be 


Ha (2 = )n ° y n 
> & 2F,( n, y—tn; 23 1-4) « a 


n=0 
Pig fatigue 8S ae 


n=0 
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For practical purposes, the terms c, are most conveniently calculated from 
the recursion formula 


2y-qt+ : 
nt = em) (Ed) Let y) ey) Hin (ey Om By —gtn— 1 xy er-1} (26) 


which is obtained from the recursion formulas for the hypergeometric functions. 
The expression (25) may serve as an example of a useful result that follows 
rather directly using the new representation of the Coulomb wave functions (11). 
It is not claimed that expansions as simple as (25) cannot be derived starting 
from standard expressions for the wave functions but in so doing it seems to 
be accomplished in a much less direct way. An attractive feature is the oc- 
currence of only the sum and difference of the linear momenta of the photons 
and electrons in the arguments and perhaps this symmetry will be helpful in 
studying more general matrix elements than the simple ones studied above. 


ACKNOWLEDGEMENT 


This work was carried out when the author had the favour of studying at the Nordic In- 
stitute for Theoretical Atomic Physics in Copenhagen. Numerical calculations of conversion 
coefficients have been performed and the author is indebted to the Nordic Institute for Theo- 
retical Atomic Physics for a journey to the computer BESK in Stockholm as well as the Swe- 
dish Atomic Committee and the BESK staff for the use of the computer. 

Accurate and fast work by Dr. 8. Hultberg facilitated the numerical work and Dr. B. 
Nagel supplied necessary and valuable criticism. 


Institute of Theoretical Physics, University of Stockholm, Stockholm, Sweden. 


REFERENCES 


1. W. Gorpon, Ann. Physik (5) 2, 1031 (1929). 

2. K. Atprer and A. WrintHer, Kgl. Danske Videnskab. Selskab. Mat. Fys. Medd. 20 nr. 18 
(1955). 

3. P. O. M. Ousson, Ark. f. Fysik. 75, 131 (1959). 

4, A, Erp&tyt, Math. Z. 40, 693 (1936). 

5. M. E. Rosr, G. H. Gorrrzer, B. I. Spryrap, J. Harr and P. Strone, Phys. Rev. 83, 79 
(1951). 

€. J. L. BurcHNatt and T. W. CHaunpy, Quart. J. of Math. (Oxford) 72, 112 (141). 

7. P. Appetit and J. Kampeé pr Fert, Fonctions Hypergéométriques et Hypersphériques. 


Gauthiers-Villars, Paris, 1926. 


Tryckt den 16 mars 1959 


Uppsala 1959. Almqvist & Wiksells Boktryckeri AB 


